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ABSTRACT 

We investigate the dynamical condensation process in a magnetized thermally bistable medium. We perform 
one-dimensional two fluid numerical simulations that describe the neutral and ionized components in the inter- 
stellar medium with purely transverse magnetic fields. We find that the clouds that are formed as a consequence 
of the thermal instability always have a magnetic field strength on the order of a few pG irrespective of the 
initial strength. This shows good agreement with the measurements of the magnetic field strength in diffuse 
clouds. We show analytically that the final magnetic field strength in the clouds is determined by the balance 
between ambipolar diffusion and the accumulation of the magnetic field due to condensation. 

Subject headings: magnetohydrodynamics — ISM: general — method: numerical — instabilities 



1. INTRODUCTION 

It is widely known that the low- and mid-temperature parts 
of the interstellar medium (ISM) are described as a thermally 
bistable medium that results from the balance of radiative 
cooling and heating (Field et al.l969; Wolfire et al. 1995, 
2003). The two stable phases are called the cold neutral 
medium (CNM), and the warm neutral medium (WNM). 

Thermal instability (TI) is the most promising formation 
mechanism of the CNM (Field 1965). Recently, many authors 
studied the dynamical condensation proc ess of the ISM driven 
by TI Koya ma & Inutsuka (2000 2002'). llnutsuka & Koyariial 
(12004 . and llnutsuka et al.l (l2005h studied the turbulent CNM 
formation as a result of TI in the layer compressed by 
shock propagation. Hennebelle & Perault (1999, 2000), 
lAudit & Hennebelld (l2005h . Heitsch et al. (2005, 2006), and 
Vazquez-Semadeni et al(2006a, b) studied an analogous pro- 
cess in supersonic conversing flows. However, there are not as 
many studies that take into account the effect of the magnetic 
field on TI (Elmegreen 1997; Hennebelle & Perault 2000; Lim 
et al. 2005; Hennebelle & Passot 2006). 

Measurements of the magnetic field in the CNM (diffuse 
clouds) from the Zeeman effect give values around 5 /iG (My- 
ers et al. 1995; Heiles & Troland 2005). It shows that thermal 
and magnetic pressures in the diffuse clouds are comparable 
(plasma beta: (3 = /^themai/ Pmag ~ !)■ Thus, the effect of mag- 
netic field on TI is obviously important. 

In the case of a one-dimensional condensation process with 
a magnetic field that is frozen in the fluid (ideal-MHD), mag- 
netic pressure increases in proportion to the square of the gas 
density, if the direction of the fluid motion is not perfectly 
aligned with the magnetic field. Thus, during condensation, 
magne tic pressure seriously affects TI. iHennebelle & Perault! 
(I2OOOI) show that magnetic pressure prevent TI, if the ini- 
tial angle between the magnetic field and the fluid velocity 
is larger than 20° - 30° with a few fiG initial field strength. 
Furthermore, even if the CNM is formed as a result of TI, 
the resulting CNM would be occupied by material with low 
f3 (see Lim et al 2005), because perfect alignment between 
the magnetic field and the velocity field is rare. This seems 
to contradict the observational suggestion that the magnetic 

' Department of Physics, Graduate Scool of Science, Kyoto University, 
Sakyo-ku, Kyoto 606-8588, Japan; tsuyoshi@tap.scphys.kyoto-u.ac.jp 
^ Department of Astronomy, University of Maryland 



field strength is independent of density in the range 0.1 to 100 
cm"-'; i.e., the WNM and the CNM have roughly the same 
field strengths (Troland & Heiles 1986). 

However, the aforementioned studies on TI that take into 
account the effect of the magnetic field do not include the re- 
alistic diffusion process of the magnetic field. In this Letter, 
to deal with the realistic process of ambipolar diffusion during 
TI, we perform one-dimensional two-fluid numerical simula- 
tions with a purely transverse magnetic field. 

2. THE MODEL 

2. 1 . Basic Equations and Numerical Scheme 

We consider the equations of one-dimensional hydrody- 
namics and magnetohydrodynamics that have purely trans- 
verse magnetic fields, including collisional coupling between 
neutral and ionized gases (Draine 1986). The neutral fluid 
part also includes radiative cooling, heating, and thermal con- 
duction. 
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where the subscripts n and / denote the neutral and ion compo- 
nents, respectively. In addition, we impose the ideal gas equa- 
tion of state pi,,i = kBp„,iT„,i/m„,i. P„ = p„ and P,- = /7, +Z?-/87r 
represent the thermal and magnetic pressure, and E„ = Pn/il- 
1 ) + p„ v^/2 and = /7,/(7 - 1 ) + p,v? /2 + B^/87r represent the 
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Table 1. Sources 



coefficient 


note 


reference 


A{T) 


3.4 X 10'5 (r/SOOOK)""' cm'g-'s-' 


1 




H^+H collisions 




K{T) 


2.5 X 10^ r" 5 erg cm-' 


2 




H+H collisions 




a{T) 


and He+ recombinations 


3 




H and He ionizations 


4 




due to CRs and X-ray 





References. — (1) Glassgold et al. 2005; (2) Parker 1953; (3) Shapiro 
: Kang 1995; (4) Wolfire et al. 1995. 



total energy. We use the specific heats ratio 7 = 5/3 and the 
mean particle masses m„ = nii = \ .21m^. The drag coefficient 
A, the recombination coefficient a, the ionization coefficient 
^, and the thermal conductivity k are chosen to simulate the 
typical ISM whose properties are listed in Table 1 . 

We use the operator-splitting technique for solving the 
equations: the second-order Godunov method (van Leer 
1979) for the inviscid fluid parts and the second-order explicit 
time integration for the cooling, heating and thermal conduc- 
tion parts. The parts of ionization, recombination, and fric- 
tional forces/heatings are updated using piecewise exact solu- 
tions for each time interval. 

2.2. Cooling Function and Numerical Settings 

We take the following simplified net cooling function which 
approximates the realistic cooling function: 

pL = n {-V + n A) erg cm"-' s~' , 
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+ 1 .4 X 10"^ \/r exp ( ^l\) 



This is obtained by fitting to the various heating (F ) and cool- 
ing (A) processes considered bv lKovama & Inutsu ka (2000). 
The resulting thermal-equilibrium state defined by C{n,T) = Q 
is shown in the bottom panel of Figure[T]as a thick solid line. 

In all of the computations shown below, we use A^ceii = 4096 
grid points covering a L = 4 pc region. Thus, the spatial res- 
olution is Ax = L/NceW — 10"^* pc, which satisfies the Field 
condition (Koyama & Inutsuka 2004). We impose a periodic 
boundary condition. The initial condition of the neutral gas is 
static and uniform thermally unstable equilibrium whose gas 
pressure is p„/k^ = 3500 K cm"-', and that of the ionized gas 
is determined from the ionization equilibrium (in this case, 

= «i/«n = 3.10 X 10~~) with the same temperature as the 
neutrals. In order to excite TI, we add density fluctuations to 
the neutral gas as 



(5pn = A ^ sin 



IttI ^ 

—x+e(i) 



(13) 



where 0{l) is a random phase. We use ^ = 0.05 Wn.ini and Zmax = 
20 in all runs. 

As shown in iHennebelle & Peraultl (l2000l) . the transverse 
component of the magnetic field strongly affects the non- 
linear evolution of TI, even if the initial angle between the 
magnetic and velocity fields is small. In order to study the 
nonlinear effect of magnetic pressure and ambipolar diffu- 
sion, we initially set a purely uniform transverse magnetic 
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Fig. 1 . — Density structure (top) and number density-pressure distribution 
(bottom) of the two-fluid case (thin solid line) and the ideal MHD case (dotted 
line) with Bi„i = 1.0 at ^ = 25.0 Myr The dashed line in the top panel 
shows the resulting ionization rate Xj = /«„ of the two-fluid calculation. 



field whose strength is Bi^i = 0, 0.03, 0.1, 0.3, and 1.0 ^G. To 
contrast the effect of ambipolar diffusion, we also perform the 
ideal MHD simulation, which has the same cooling function 
and thermal conduction with Bi„i =1.0 fiG. 

3. RESULTS 
3.1. The density evolution 

In all runs except the ideal MHD case, TI can fully develop 
and can create a two-phase medium (CNM and WNM). On 
the other hand, in the ideal MHD case, the collapse of the 
dense region stops before the density reaches that of the ther- 
mally stable CNM branch because the magnetic pressure be- 
comes strong enough to support the cold region. Figure [T] 
shows the resulting density structure (top panel) and number 
density-pressure distribution (bottom panel) of the two-fluid 
case (thin solid line) and the ideal MHD case (dotted line) 
with Bini = 1 .0 /iG at f = 25.0 Myr 

In the case of two-fluid calculation, de nsity evolution of the 
neutral gas is similar to that reported in iKovama & Inutsuka 
(I2OO6I) . At first over-density perturbations grow to the CNM 
phase. The most unstable scale of the perturbation is writ- 
ten as Iji ~ ^/lJF (Field 1965) that is ~ 1 pc in the unsta- 
ble equilibrium and 0.1 pc in the CNM, where /(, = Csfcooi is 
a acoustic scale (the sound traveling scale within a cooling 
time) and /p = \/ kT / n^A is the Field length. The growth 
timescale of the condensation is determined by the cooling 
time fgrow ^ fcooi = T /nA/{j- 1) ~ 1 Myr in the unsta- 
ble equilibrium and ~ 0.1 Myr in the CNM. Thus the growth 
timescale becomes faster as they condense. Figure |2] shows 
the neutral gas density evolutions of the densest point in the 
range x G [0.0, 1.5] pc in which contains only one evolving 
density peak for all run. Clearly in the all two-fluid case the 
TI growth as fast as the unmagnetized case {Bi„i = 0). This 
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Fig. 2. — Neutral gas density evolutions of the densest point in the range 
a:G [0.0, 1.5] pc. 
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Fig. 3. — Magnetic field strength evolutions of the densest point in the 
ranges: G [0.0, 1.5] pc. 



shows that the ambipolar diffusion effectively drift magnetic 
field from the condensing region. A quantitative estimation of 
a critical magnetic field strength that the ambipolar diffusion 
becomes effective is shown in the next section. 

After the growth of over-density perturbation, less-dense 
perturbations grow to the WNM phase, because growth 
timescale of the rarefaction is determined by the heating time 
that is longer than the cooling time (fheat ~ 1 Myr in the un- 
stable equilibrium and ^ 10 Myr in the WNM). 

3.2. The pressure evolution 

During the evolution, the total pressure pxai = Pn + Pi + 
B^/Stt is spatially uniform owing to the fact that the most un- 
stable scale of the TI is smaller than the acoustic scale /ji < L- 
The condensation terminates when the neutral gas reach the 
stable phase or magnetic pressure becomes larger than the 
neutral gas pressure. In our calculations, all two-fluid runs 
corresponds to the former case, and the ideal-MHD run cor- 
responds to the later case. 

3.3. The magnetic field strength evolution 

Figure[3]shows the magnetic field strength evolutions of the 
densest point in the range x e [0.0, 1 .5] pc. Amazingly, mag- 
netic field strength in the CNM always becomes on the order 
of a few fiG irrespective of the initial strength! The resulting 
p in that points at f = 10 Myr are /3 = 2.46, 5.55, 12.2, and 
28.2 that correspond to the runs Bini = 1.0, 0.3, 0.1, and 0.03 
/iG, respectively. As shown in the next section, this results 
can be understood that the critical strength that the ambipolar 
diffusion sets in becomes on the order of a few fiG, when the 
CNM is generated by the TI. 

4. CRITICAL MAGNETIC FIELD STRENGTH 

We can estimate the condensation speed of the TI as 



the drift speed can be written as 



V'TI ■ 



^cool 



(14) 



where the spatial- and time-scale are essentially determined 
by the most unstable scale of the TI and the cooling time, 
respectively. 

The drift speed of the magnetic fieldline from condensing 
region can be estimated as follows. During the condensa- 
tion the ionization degree is small {Xi = 10"^- 10"^), thus the 
friction- and the Lorentz-forces are in balance. In that case 



V'drift ^ 



|(V xE)xB\ 
47rApnPi ^ 47r AXipl hi' 



(15) 



(Nakano 1984; Shu 1992) where we simply replace the spatial 
derivative with the inverse of the most unstable scale of the TI 
that is the characteristic scale of the system. 

If the drift speed is slower than the condensation speed 
^drift < v-Yi, magnetic fieldlines are also condensed with neu- 
trals. Such a case realizes when the magnetic field is weak. If 
the drift speed is faster than the condensation speed v'diift > vti, 
magnetic fieldlines drift from condensing region. Such a case 
realizes when the magnetic field is strong. Therefore, during 
the evolution, the magnetic field strength inside the condens- 
ing region is adjusted to Vdritt ^ vji- 

Equating equations ( fT4l ) and ( fTSl l. the critical strength that 
the ambipolar diffusion becomes effective is written as fol- 
lows 

BAD=^^TrAXiplcjF. (16) 

If we assume the ionization equilibrium and isobaric evolu- 
tion, equation ( fT6] l is written as a function of the neutral num- 
ber density that is plotted in Figure |4] Note that the ioniza- 
tion equilibrium is not maintained during the condensation 
because the timescales of the recombination and ionization 
are, respectively, slightly smaller and larger than the cooling 
time. Thus this estimation overestimates the critical value dur- 
ing condensation. This figure shows that the critical strength 
Bad is not sensitive to the density. Therefore we can express 
the value of the critical strength using typical values in the 
CNM: 



Ba 



:3.2 



Pn/ks 



4000 K cm 



1.475 



50 cm- 



-0.4875 



/iG, (17) 



where we use ^ ~ 1.0 x 10 
-^andA~3.6x 10-2'*r 



-1, a -6.8 X 10-11 cm^ 
and A ~ i.b x W'-' T" " erg cm^^ (Wolfire et al. 2005). 
This value shows good agreement with our numerical experi- 
ments and the measurements in the diffuse HI cloud, and that 
is why the resulting magnetic field strength in the CNM al- 
ways converge to the order of a few /iG irrespective of the 
initial strength. 

5. SUMMARY AND DISCUSSIONS 

We studied dynamical condensation process of the TI that 
takes into account the effect of magnetic pressure and the 
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Fig. 4. — Critical magnetic field strength vs. number density of the neu- 
trals. 



ambipolar diffusion by one-dimensional two-fluid numerical 
simulation. We showed that even if /3 ~ 1 initially, the TI can 
fully develop and generate the moderately magnetized cloud 
(/? ~ 1) owing to the ambipolar diffusion, whereas prevented 
from condensation due to magnetic pressure in the case of the 
ideal-MHD. Therefore, if diffuse HI clouds are formed by the 
TI, it is natural that observations show no increase of magnetic 
field strengths over the density range 0.1 to 100 cm~''(Troland 
&Heiles 1986). 



Furthermore we showed that magnetic field strength in the 
CNM converges to the order of a few /iG irrespective of the 
initial strength. This result can be understood that during the 
condensation the magnetic field strength is tuned to the criti- 
cal value owing to the accumulation of magnetic fieldline by 
the TI and ambipolar diffusion. The critical value becomes 
around 3 /iG in the typical ISM (see equation ITtI ). This 
value shows good agreement with the measurements of the 
magnetic field strength in the diffuse HI clouds (Myers et al. 
1995; Heiles & Troland 2005). Thus our result suggest that 
even if there are large fluctuation of magnetic field strength in 
the WNM, the strength in the diffuse clouds always becomes 
a few fiG. 

In this study, for simplisity, we choosed initial condition to 
the thermally unstable equilibrium. A study of more realistic 
situation, for example converging two flows or shock propa- 
gation, will be explored in our next paper. In the converging 
two flows case, ram pressure of the flows prevent ions drift, so 
ambipolar diffusion can become effective after the converging 
flows become weak. 
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